Magnetic Phase Diagram of Ca2_xSra;Ru04 Governed by Structural Distortions 
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We constructed, by the first-principles calculations, a magnetic phase diagram of Sr2Ru04 in the 
space spanned by structural distortions. Our phase diagram can qualitatively explain the experi- 
mental one for Ca2-a;Sra;Ru04. We found that the rotation and the tilting of RuOg octahedron are 
responsible for the ferro- and antiferro-magnetism, respectively, while the flattening of RuOe is the 
key factor to stabilize those magnetic ground states. Our results imply that the magnetic and the 
structural instabilities in Sr2Ru04 are closely correlated cooperatively rather than competitively. 
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FAGS numbers; 75.30.Kz, 74.70.-b, 71.27.+a 

Both the magnetic and the structural instabilities are 
essential issues for the unconventional superconductivity 
in Sr2Ru04 which is the only example of a non- 

cuprate layered perovskite superconductor. It was first 
suggested that the Sr2Ru04 is close to the ferromagnetic 
(FM) instability with strong FM spin fluctuations, 
which may naturally lead to a spin-triplet p-wave pairing 
mechanism ]7|-[To|]. However, the recent observation | |Tl| ] 
of sizable antiferromagnetic (AF) incommensurate spin 
fluctuation, due to the Fermi surface nesting in- 
dicates that more careful studies are needed. As for 
the structural aspect, it was pointed out by experiments 
that Sr2Ru04 is very close to the structural instability 
with respect to the RuOg rotation With such a 

situation, one may consider that three kinds of insta- 
bilities, superconducting, magnetic and structural ones, 
may compete. Nevertheless, the correlation among those 
instabilities has not been fully discussed. It was found 
recently |Q that the cleaved surface of this material is 
reconstructed to form the c(2x2) structure which can be 
regarded as the frozen RuOg rotation mentioned above. 
Furthermore, the density-functional calculation predicts 
that the surface ferromagnctism is strongly stabilized by 
the structural reconstruction [Q. This prediction sug- 
gests that the structural and magnetic instabilities co- 
operate rather than compete, although the surface ferro- 
magnetism has not been experimentally confirmed up to 
now. 

On the other hand, the recent studies on the 
Ca2-2;Sr2;Ru04 suggest strongly the cooperative feature 
of the structural and magnetic instabilities in the bulk. 
Moreover, the system shows a very rich phase diagram 
and provides us with an opportunity to analyze the cor- 
relation between the magnetism and the structure more 
extensively. Below is a brief description of the experi- 
mental observation for Ca2-2;Sr2;Ru04 by Nakatsuji et 
al. With the Ca substitution for Sr, the system 

is successively driven from the non-magnetic (NM) 2- 
dimensional (2D) Fermi liquid [x ^ 2.0) to a nearly 
FM metal {x ~ 0.5), an antiferromagnetically corre- 



lated metal (0.2 < a; < 0.5), and finally an AF insulator 
(x < 0.2). Since the substitution is isovalent, the dom- 
inant effects are the structural modifications due to the 
reduced ionic size of Ca compared with Sr. Evidence |l^ ] 
has been presented by neutron scattering that the struc- 
tural distortions caused by the Ca substitution correlate 
with the changes in the magnetic and the electronic prop- 
erties. 

The main aim of this letter is to study how and why 
the magnetism of Ca2-a:Sra:Ru04 is affected by struc- 
tural distortions. In order to extract essential aspects, 
we assume that for a given crystal structure, the elec- 
tronic structure is not affected by the relative content 
of Ca and Sr. Therefore, in the following, we study 
the stable magnetic phases of Sr2Ru04 for given struc- 
tural distortions. Three types of structure distortions, i.e. 
RuOg octahedron rotation about the c-axis, RuOg tilting 
around an axis parallel to the edge of octahedron basal 
plane and the flattening of RuOg along the c-axis are 
identified from experiments p^ ] . Our phase diagram can 
qualitatively explain the experimental phase diagram of 
Ca2-2;Sr3;Ru04, demonstrating the crucial roles of struc- 
tural distortions for the tuning of electronic and magnetic 
properties, and further supporting our previous predic- 
tion for the surface . In particular, we found that the 
RuOg rotation can enhance the FM instability signifi- 
cantly, while the combination of tilting and rotation of 
RuOg is responsible for the enhancement of AF insta- 
bility. Furthermore, we point out that the fiattening of 
RuOg is a key factor to stabilize the magnetic (both FM 
and AF) ground states. The basic physics governing the 
phase diagram can be understood in terms of the strong 
coupling between the lattice and the magnetism through 
the orbital degrees of freedom. Our results strongly sug- 
gest that, in Sr2Ru04, the magnetic fluctuations can be 
significantly enhanced by the structural fluctuations, im- 
plying the necessity of reconsidering the coupling mech- 
anism in the bulk superconductivity. 

The calculations were performed with the first- 
principles plane-wave basis pseudopotential method 
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based on the local density approximation (LDA). The va- 
lidity of LDA treatment for ruthenates was demonstrated 
in Ref. The 2p states of oxygen and 4d states of 

Ru are treated by the Vanderbilt ultrasoft pseudopoten- 
tial while the norm-conserving scheme is used 
for other states. The cutoff energy for the wave function 
expansion is 30 Ry. The A:-point sampling of the Brillouin 
zone (BZ) was well checked to provide enough precision 
in the calculated total energies. The theoretically opti- 
mized lattice parameters a = 3.84 A and c — 12.70 A for 
the bulk Sr2Ru04 are in good agreement with the exper- 
imental ones a = 3.86 A and c 12.73 A. The degree of 
flattening of RuOg octahedron A is defined by A = dc/dab 
with dc [dab) denoting the Ru-0 bond length along the 
c-axis (in the afc-plane) with the RuOg volume fixed. 
Rotation and tilting of RuOg octahedron are operated 
with the Ru-0 bond lengths fixed. In order to construct 
the magnetic phase diagram, the lowest energy magnetic 
phase for each crystal structure is searched for among 
different (NM, FM and AF) phases. In the present work, 
we focus our attention only on phases described within 
the c(2x2) unit cell. The soft phonon mode of S3 at 
the zone boundary in Sr2Ru04 ||l^ and the AF state of 
Ca2Ru04 are in this category. 

Figure 1 shows the calculated phase diagram |2^ , while 
the Table 1 summarizes the calculated total energies 
and magnetic moments for some particular points in the 
phase diagram. Hereafter, <j) and 6 denote the rotation 
angle and the tilting angle, respectively. The apical oxy- 
gen and the oxygen in the a6-plane are called 0(2) and 
0(1) respectively. From right to left of the phase dia- 
gram, first the RuOg starts to rotate along the c-axis by 
up to 12°, and then with the 12° rotation being fixed, the 
RuOg starts to tilt up to 12°. The structural analysis by 
the neutron scattering jl^ allows us to make a one-to- 
one correspondence between the structural changes, i.e. 
the horizontal axis of our phase diagram, and the doping 
level X in Ca2_2;Sr2;Ru04. For x = 2.0 (Sr2Ru04), the 
system has lA/mmm symmetry with (j) ^ 9 = 0°, corre- 
sponding to the right end of our phase diagram. With 
reduction of x, RuOg starts to rotate and the symmetry 
is reduced to Mi/acd until x — 0.5 (Cai.5Sro.5Ru04), 
where = 12.78° and 6* = 0° at 10 K. With further re- 
duction of x, RuOg starts to tilt and the symmetry is 
further reduced to Pbca until x — 0.0 (Ca2Ru04) where 
(j) = 11.93° and 9 ~ 12° at low temperature, correspond- 
ing to the left end of our phase diagram. It was also 
pointed out by the experiment p6| that, from x = 2.0 
to X = 0.5, the degree of flattening A remains almost 
constant (~ 1-07), while from x — 0.5 to a; = 0.0, the ro- 
tation angle (p is almost unchanged 12°). Three rep- 
resentative experimental points are shown in our phase 
diagram by triangles. Now, the basic tendency suggested 
by our phase diagram is that the RuOg rotation will drive 
the system from a NM state to a FM state, while the 
subsequent tilting plus the flattening of RuOe will push 



the system to an AF region. This general tendency is 
quite consistent with the experimental results. It should 
be noted that the rich experimental phase diagram can 
be simply understood in terms of the close coupling be- 
tween structural distortions and magnetism. Although 
the real long range FM ordering in Cai.5Sro.5Ru04 is not 
confirmed yet, the significant enhancement of spin sus- 
ceptibility in this doping level undoubtedly implies the 
strengthening of FM instability. Another important as- 
pect in our phase diagram is that the flattening of RuOg 
is so important not only for the AF state but also for the 
FM state. This suggests that simply by uniaxial pressure, 
the Sr2RuOi can be driven from the NM state to a FM 
state. 

The basic questions concerning our phase diagram are: 
1) Why are the RuOg rotation and tilting correlated with 
the tendency to the FM and AF states? 2) Why is the 
RuOg flattening so important for the magnetic solutions? 
Before answering these questions, let us discuss the role 
of each Ad orbital in the electronic properties of Sr2Ru04, 
which is essential to our later discussions. The three Ru 
t2g orbitals (dxy, dyz^ dzx) hybridize with each other only 
very weakly in tetragonal Sr2Ru04. Therefore each or- 
bital plays distinct roles. The projected density of states 
(DOS) shown in Fig. 2(a) indicates that the d^y orbital 
contributes dominantly to the well-known van Hove sin- 
gularity (VHS) just above the Fermi level. The 7 Fermi 
surface has the character of dxy It is mostly responsible 
to FM spin fluctuation due to the high DOS around the 
Fermi level. On the other hand, dyz and dzx orbitals con- 
tribute to the a and [3 Fermi surfaces and produce the 
incommensurate spin fluctuation coming from the strong 
nesting effect due to the quasi-one-dimensional nature 
of those states. The calculated bare spin susceptibility 
shown in Fig. 3(a) for the undistorted compound has 
the incommensurate peak at Q=(27r/3a, 27r/3a), being 
consistent with previous calculations [jl^ . In real mate- 
rials, the two factors, i.e. the FM instability due to the 
high DOS at the Fermi level combined with the q depen- 
dent Stoner factor J@,|2| and the AF instability due to 
the nesting effect [ [l2[ , will compete. 

The RuOg rotation couples mostly with the dxy orbital 
but not with the dyz, dzx orbitals because the pdir type 
hybridization between the 0(l)-2p and the dxy states will 
be significantly reduced by the RuOe rotation, but those 
between the 0-2p and the dyz , dzx states are not affected 
so much. The direct results of this reduced pdn type 
hybridization between the 0(l)-2p and the dxy states 
arc, first the narrowing of dxy band width and second 
the downward shift of dxy band, as shown in Fig. 2(b) 
(about 0.4 eV narrowing and 0.1 eV downward shift of 
dxy band for (f) = 12°, A — 1.07). As the latter brings the 
VHS closer to the Fermi level, both of the two results 
will enhance the DOS at the Fermi level. Another effect 
coming from the downward shift of the dxy states is the 
population reduction in the dyz, dzx states, which may 
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shift the Fermi surface nesting vector closer to the zone 
boundary. However, the increase of the DOS at the Fermi 
level is the dominant effect and the tendency towards FM 
instability is enhanced by the RuOg rotation. Once tilt- 
ing is additionally introduced, all of the t2g bands will 
become narrower. This will enhance the nesting effect 
and enhance the AF instability. The discussion so far 
can answer the first question. 

Now let us discuss the effects of RuOg flattening. 
There are two factors also. First, with the flattening 
of RuOe octahedron, the Ru-O(l) bond length will in- 
crease, while the Ru-0(1)-Ru angle remains 180°. The 
increased bond length will reduce all the pdir type hy- 
bridizations between the 0(l)-2p and the d^y, dyz, d^x 
states. Therefore, width of all the bands of three Ru-t2g 
states is reduced as shown in Fig. 2(c) (about 0.4 eV for 
dyz,zx bands and 0.3 eV for d^y band for A = 0.96), mak- 
ing the DOS at the Fermi level higher. This will favor 
the FM solution. Another very important results of flat- 
tening is the orbital polarization. It is obvious that the 
tetragonal distortion by the flattening will populate the 
dxy states and depopulate the dyz, dzx states (about 0.2 
eV downward shift of dxy band for A = 0.96). The effect 
is similar to the RuOg rotation. The orbital polarization 
due to flattening will also shift the nesting vector to the 
zone boundary as shown in the susceptibility calculations 
(Fig. 3). This will favor the commensurate AF state of the 
system. The net effect by rotation, tilting and flatten- 
ing of the RuOg will depend on the competition among 
them and the phase diagram of Fig.l demonstrates the 
situation in a space spanned by those distortion modes. 
Figure 2(f) shows the calculated DOS for the AF state 
with (j) ^ 9 — 12°, i.e. almost the experimental structure 
of Ca2Ru04. It is clear in this case that the occupied mi- 
nority spin states mostly come from the dxy orbital due to 
flattening of RuOg. Therefore, the strong superexchange 
interaction between the occupied majority-spin and un- 
occupied minority-spin dyz, dzx orbitals will stabilize the 
AF ground state 

In summary, by constructing a phase diagram of 
Sr2Ru04 with structural distortions, we find the strong 
coupling between the lattice and the magnetism. Our 
phase diagram can qualitatively explain the experimen- 
tal phase diagram of Ca2-xSr2;Ru04. We demonstrate 
that the RuOe rotation will enhance the FM instabil- 
ity in the system, while the tilting plus the flattening of 
RuOe make the system AF. We pointed out that the flat- 
tening of RuOg is so important not only for the AF state 
but also for the FM state. An important implication of 
our results is that the magnetic and the structural insta- 
bilities in Sr2Ru04 should be strongly correlated. The 
structure fluctuation and the magnetic fluctuation coop- 
erate. Actually the phonon mode corresponding to the 
RuOg rotation is quite soft |]l^ in the bulk, and this ro- 
tation will enhance the FM instability. All these results 
imply the necessity of reconsidering the coupling mech- 



anism for the unconventional superconductivity. In this 
context, we propose a possible way to identify experimen- 
tally the relationship between the FM fluctuation and 
the superconducting state. As the uniaxial compression 
of Sr2Ru04 will enhance the FM fluctuation without in- 
troducing the disorder, the variation of superconducting 
transition temperature against uniaxial compression may 
provide important information. 
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FIG. 3. A contour plot of the calculated bare spin suscep- 
tibility for {a)cj)=0=O°, X = 1.07; {h)(f)=0=O°, X = 0.96. The 
dashed lines are guide to the eye for the nesting vectors. The 
red color denote the higher intensity. 
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TABLE I. The calculated total energies and magnetic mo- 
ments for some particular points in our phase diagram. 





NM 


FM 


AF 


A = 0.96 


meV 


-40 meV 


-8 meV 


(t> = e = Q° 




1.26 mb/Ru 


0.7 i^b/Ru 


A = 1.07 


meV 


-25 meV 




0= 12°, 6' = 0° 




0.74 ixb/Ru 




A = 0.96 


meV 


-100 meV 


-117 meV 


^ = 12°, e = 12° 




1.13 /ub/Ru 


0.93 //s/Ru 



